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ABSTRACT 

We model the latest HST WFPC3/IR observations of > 100 galaxies at redshifts z = 7-8 in terms 
of a hierarchical galaxy formation model with starburst activity. Our model provides a distribution 
of UV luminosities per dark matter halo of a given mass and a natural explanation for the fraction 
of halos hosting galaxies. The observed luminosity function is best fit with a minimum halo mass per 

galaxy of 10 -° 9 M Q , corresponding to a virial temperature of 10 -°- 7 K. Extrapolating to faint, 
undetected galaxies, the total production rate of ionizing radiation depends critically on this minimum 
mass. Future measurements with JWST should determine whether the entire galaxy population can 
comfortably account for the UV background required to keep the intergalactic medium ionized. 

Subject headings: Cosmology: theory — early universe — galaxies: formation — galaxies: high-redshift 
— stars: formation 



1. INTRODUCTION 

The pursuit for the first galaxies has recently entered 
a new phase as observations at redshifts z > 6 have now 
probed the epoch before cosmic reionization was com- 
plete enough to fill in the Gunn-Peterson trough. The 
WFC3/IR Camera on HST has recently detected a sam- 
ple of more than 60 z ~ 7 and nearly 50 z ~ 8 Lyman- 
break galaxies (LBGs) (|Bouwens etai]|2010d) and pro- 
vided c onstraints on the gala xy abundance as early as 
z ~ 10 (jBouwens et al.ll2010b| ). The UV spectral slopes 
of these faint sources were found to be very flat, per- 
haps i ndicating dust-free environments (|Bouwens et al.l 
l2010al) . and their stellar masses have be en inferred from 
measurements in the rest -frame optical (Gonzale z et al.l 
20101: lLabbe"et~aT]l2010al |b1). These observations inform 
theoretical models of galaxy formation and attempt to 
probe the amount of radiation available to affect the ion- 
ization state of the intergalactic medium (IGM) but are 
often limited by survey sensitivity. The traditional in- 
terpretation of such data relies on an assumed ratio of 
UV luminosity to star formation rate (SFR) that requires 
burst ages longer than the exponen tial burst time-scal e 
and time-scales longer than 1 Gyr {Madau et "all [l998h . 
This assumption cannot be satisfied at redshifts z > 6, 
where the age of the Universe is shorter than 1 Gyr. 

Many theoretical studies based on numerical and semi- 
analytic techniques, have shown that the luminosity 
function (LF) of high-redshift LBGs can be explained 
by the hierarchical formation of dark matter halos 
whose associated bar yonic gas forms stars in merger- 
generated bursts (e.g. [Baugh et al.l |2005| iFinlator et al.1 
I2010t lLacev et all 120101 : iSalvaterra et al l 12010ft AnaT 
lytic work has tried to fit simpler models to the ob- 
served LF to probe the d uty cycle and mass-to-light ra- 
tio o f observed galaxies (jStark et al.l 2007; Trenti et al.1 
120101 e.g) assigning a single galaxy luminosity for each 
halo mass, while others have considered that the mass 
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of a host halo may merely define the probability dis- 
tribut ion from which a galaxy's lu m inosity is drawn 
(e.g. ICoorav fe Milosavlievid l2005a|Tbl: iCoorav fc Ouchil 
20061: iVale fc Ostrikerl 12004 120061 120081) Associating 
the mass of underlying halos with observed luminosity is 
crucial for describing the clustering properties and bias 
of high-redshift galaxies as well as the contribution of 
cosmic variance to fluctuations in the measured abun- 
dance from field-to-f i eld (e.g. iTrenti fe Stiavelll 120081: 
Munoz fc Loebll2008at lOverzier et al.ll2009t iMunoz et alii 
2010t iRobertsonl l2010al lbTr The relationship may also 
provide insights into how much ionizing radiation is pro- 
vided by galaxies too faint to be detected with current 
instruments. 

The amount of currently undetected UV radiation at 
high redshifts is unknown. While the Early Release Sci- 
ence observati ons with WFC3/IR c an probe down to 
27.5 AB mag (Bouwcn s et al.l [2010d ). it is almost cer- 
tain that man y fainter galaxies remain to be observed 
bv J WST (e.g. IBarkana fc Loebl l2000bl: IWvithe fc Loebl 
120061 : Salvatcr ra fc Ferrara!l2006| ). Galaxies should exist 
in halos down to a mass below which the assembly or re- 
tainment of gas is suppressed. While these dwarf galaxies 
near the suppression limit may be faint, their abundance 
may make them, in aggregate, large contributors to the 
total UV background. The exact suppression mass of 
galaxy formation is due to an unknown combination of 
the heating of the IGM during reionization and thermal 
and mechan ical feedback by intern al mechanisms such as 
supernovae (IWvithe fc Loebll200l . 

We propose that the suppression of galaxy formation 
below a certain halo mass threshold may already be evi- 
dent from the currently observed LFs at z > 6. Because 
bright galaxies are formed primarily through mergers 
among fainter ones and since the contribution to the pop- 
ulation of bright galaxies from lower mass halos shining 
at the bright end of their luminosity distributions need 
not be negligible, we expect the decrease in the faint end 
of the LF due to the suppression of galaxy formation to 
be gradual and extend to larger luminosities than previ- 
ously anticipated. We couple hierarchical merger trees to 
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a simple model of star formation to calculate the luminos- 
ity distribution function (LDF) for galaxies as a function 
of their host halo mass and the resulting LF. The sup- 
pression mass of galaxy formation applied to the variety 
of merger histories provides us with a physically moti- 
vated explanation for the fraction of halos that are form- 
ing stars at the time of observation (i.e. the duty cycle) 
as a function of mass. We calculate the amount of unob- 
served UV rad iation at z = 6 , 7, an d 8 and test the appli- 
cability of the Ma dau et al.l (|1998f ) relati onship between 
UV lu minosity and SFR. We note that iRaicevic et al.l 
(2010) recently considered the effect of photoionization 
on t he LF using the GALFORM model of galaxy forma- 
tion (|Cole et alJl2000l:IBaugh et~aT1[2005l) . but our models 
produce very different results at the faintest luminosities. 
This is likely because photoionization in the GALFORM 
model only affects gas cooling while allowing already cold 
gas to form stars even in the smallest halos. 

In ^ we describe the merger trees, star formation 
model, and suppression prescriptions used to populate 
the LDF for each halo mass. We describe the resulting 
shape and mass-dependence of the LDF in $3] and dis- 
cuss the physical origin of the luminous duty cycle of 
halos in S|4] I n $SJ we fit the mean of the LDFs, the re- 
sulting star formation efficiency, and the mass at which 
star formation is suppressed to match the observed LFs 
at z > 6. In <J6] we compare the star formation rate to 
the prediction from the UV luminosity. We then dis- 
cuss the implications for the ionization state of the IGM 
of the suppressed star formation in low-mass halos and 
the abundance of faint galaxies yet to be observed in SJ71 
Finally, we summarize our main conclusions in <J5] 

2. THE MODEL 

Our calculation of the LDF has two main compo- 
nents: a merger tree builder and a star formation model. 
Throughout, we assume a flat, ACDM cosmology with 
cosmolo gical parameters taken from the WMAP- 5 data 
release (jKomatsu et al.ll2009t iDunklev et al.ll2009t l. 

2.1. Merger Histories 

We generate merger trees based o n the extended Press- 
Schech ter procedure outlined in iSomerville &: Kolattl 
(1999). The method selects, for each descendent, a series 
of progenitors from the mass- weighted conditional mass 
function truncated at an upper mass limit for each sub- 
sequently selected progenitor such that the total mass in 
progenitors does not exceed the descendent mass. Once 
the difference in mass between the descendent and the 
growing list of progenitors falls below M res , the resolution 
limit of the algorithm, the remaining mass is assigned as 
diffuse accretion. If a descendent has two or more pro- 
genitors, we determine the merger ratio by considering 
the two largest progenitors. If the mass ratio between the 
two largest progenitors is smaller than 1:3, we denote the 
interaction to be a major merger; all other configurations 
are minor mergers. We have tested our procedure using a 
threshold ratio of 1:7 and found no noticeable difference 
in our results. The algorithm is then iterated on each 
progenitor until the tree is ended when the masses of all 
progenitors have fallen below M res . 

2.2. Starbursts 



Each halo in the merger tree is assumed initially to con- 
tain an amount of baryonic gas equal to a fixed fraction 
(f2b/f2 m ) of its halo mass. This gas is gradually converted 
into stars through bursts of star formation. There is a 
great body of evidence that the starbursts that illuminate 
LBGs are generated in mergers rather than in a q uiescent 
mode fe.g. lBaugh et alJl2005HLacev et al.ll2010D . There- 
fore, we ignore all quiescent star formation and generate 
starbursts exclusively during major mergers. 

After a minor merger, all starbursts taking place in 
associated branches of the tree are allowed to continue 
simultaneously with their own reservoirs of gas. If not 
fully coalesced at the time of observation, these simul- 
taneous bursts may appears as multiple cores in the 
galaxy morpho logy or simply be b eyond our current abil- 
ity to resolve (jOesch et al.l I2010T ) . However, in a major 
merger, all of these bursts are shut off and a new one 
is begun. The gas remaining from all progenitors is as- 
sumed to be instantaneously funne led to the center w here 
it forms a new d isk. F ollowing IMo et al.l (jl998l) and 
iBarkana fc Loebl (|2000aD . we assume the disk to have 
an exponential shape such that the surface density falls 
At high redshift, when the energy 
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density of the universe is dominated by the contribution 
from matter, the corresponding exponential size scale of 
the disk is given by 
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where we take the specific angular momentum of the disk 
to be equal to that of the halo (i.e. jd/ m d = 1), r vir 
is the halo virial radius, v c is the circular velocity of 
the halo, H is the Hubble parameter, and A is the spin 
parameter which we draw randomly for each disk from 
a log-normal distribution centered at A = 0.05 with a 
standard deviation a\ = 0.5 in log-space. The central 
surface den sity is E n = M £a j 2 ir R^. 

Following iKennicuttl (|1998l) . the surface star formation 
rate density is given by Ssfr = eSgasAdyni wnere e is 
the fraction of gas converted to stars in a dynamical 
time idyn- We verified that the disks under considera- 
tion are unstable to fragmention (i.e. have a Toomre 
Q-parameter smaller than unity). Using surface densi- 
ties averaged over the exponential scale radius of the disk 
and idyn = 2 ir R^/v c , e is found empirically to be ~ 20%. 
However, this relation also holds in azimuthally-averaged 
rings at radius r with td yn = 2itr/v c . Since we are in- 
terested in the total star formation rate produced entire 
disk from gas added at any radius, we integrate through 
the disk considering separately the contributions inside 
and outside the exponential scale radius. 
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Substituting for i?d and So, we find 
M*«0.66M Q /yr 
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Since M* = — M gas in a single burst between major merg- 
ers, equation ([3]) represents a differential equation in M* 
(or M gas ) whose solution is an exponential with a time 
scale given by: 
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where 



1 + z 
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(4) 



(5) 



is the age of the universe at redshift z. 

In each time step, we use results fro m a simple stel- 
lar po pulation generated by Starburst99 (|Leitherer et al.1 
1999) to calculate the contribution of each newly added 
group of stars to the final luminosity at 1500 A. We as- 
sume a Salpeter initial mass function (IMF) with a slope 
of 2.35 between 1 and 100 M Q and a metallicity 4% of 
solar. 

Finally, we include the effect of a suppression in galaxy 
formation below halos of a given mass. Some combina- 
tion of processes, such as supernovae feedback or pho- 
toionization, may push or heat the gas so that it escapes 
the gravitational potential well of the halo. Thus, no 
starbursts will be generated in halos smaller than M supp , 
and neither can a halo smaller than M sup p be a con- 
stituent in a starburst-generating major merger even if 
the resulting halo is larger than M supp . These two condi- 
tions (but especially the latter) combine to make inactive 
even some halos larger than M supp (see Q ■ The simplest 
way to incorporate these effects into our model is simply 



to set M r , 



Msi 



in generating the merger trees. This 



prescription is not quite realistic, since the feedback pro- 
cesses that suppress star formation are undoubtedly time 
dependent, especially during reionization. However, for 
simplicity, we assume that the contributions to the LF 
from minihalos smaller than M supp before reionization 
are minimal by the redshifts considered here. For conve- 
nience, we define 



m h = logio (M h /M ), 
m supp = logio (M supp /M ). 

3. THE LUMINOSITY DISTRIBUTION FUNCTION 



(6) 



Our model results in an approximately log-normal dis- 
tribution for the UV galaxy luminosities (1500 A) pro- 
duced by halos of a given mass: 
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in agreement with p revious as sumptions (e.g. 
Cooray fc Milosavlievicl |2005a|li]; iCoorav fc Ouchil 



2006tlVale fc Ostrikerll20o4 12006020081 ) . As anticipated 
by th e self-similarity of halo mergers (|Fakhouri et al.l 
l2010f ). we find that, independent of redshift, L 



proportional to halo mass. iBouwens et~aH (|2008f ) 
pr eviously estimated a power-law slope of 1.24 based 
on IBouwens et al.l (|2007l ) data at z ~ 4. We reiterate, 
that many of the assumptions that went into our 
model, including the lack of a quiescent component 
to star formation and dust extinction, are only valid 
at redshifts beyond 6. We do not find a change in 
the proportionality of lu minosity to halo mass at high 
masses as considered by IBouwens et al.l (|2008f) . Since 
the timescale for the coalescence of subhalo s after merger 
is re l ated to the ratio of their masses (W etzel et al.l 
[200llWetzel fc WhitdlMlCh and we have selected major 
merges based on a fixed mass ratio, we do not expect 
a fall off in the rate of major mergers in more massive 
halos. Finally, we also find <7l between 0.2 and 0.3 
for all masses and redshifts considered. Consequently, 
the galactic luminosity produced by halos of the same 
mass can easily vary by ~ 1.5 magnitudes or more. For 
convenience, we will assume a fixed value of ctl = 0.25 
for all calculations of the LF throughout the rest of this 
paper. 

4. THE LUMINOUS DUTY CYCLE OF HALOS 

We first clarify a slight ambiguity of definition in the 
literature. The duty cycle, enc, is defined as the fraction 
of a halo's lifetime over which it is luminous. If halos fluc- 
tuate stochastically on and off, then cdc also represents 
the probability that a halo will be on at the moment of 
observation and the fraction of all halos at that time that 
are on. Considering halos of a given mass for whom a 
single luminosity has been assumed, this concept of the 
duty cycle lowers the abundance of halos observed from 
that predicted in the halo mass function, i.e. 



dn ohs (L) 
dL 



dn dL(M h ) 



dM h dM h 
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is 



In various models, eoc m ay be a function of variables 
such as mass or redshift or left as a constant parameter. 

However, if the luminosity of a halo or the frequency 
of its being in a luminous state is not constant in time, 
the fraction of observable halos for a given halo mass 
need not be equal to cdc- I n our model, there are two 
reasons why a halo may not be observable. The first 
is that, given the continuous distribution in luminosity 
for a given halo mass, some are not bright enough to be 
detectable at their distance. However, these halos will 
simply appear in another luminosity bin, and we will 
proceed by first calculating the full LF and subsequently 
applying an observable threshold for a given survey. The 
second reason why a halo may not be observable is be- 
cause, in the limited history of the universe at the mo- 
ment of observation, its merger tree does not contain a 
single major merger whose constituents were more mas- 
sive than M supp . Thus, according to our model, it will 
not have had even one starburst, and will remain com- 
pletely dark. Halos much more massive than M supp have 
had at least one starburst-generating merger, while those 
closer in mass to M supp may not have since many of its 
recent progenitors are below M supp . We define the prob- 
ability that a halo of a given mass has had at least one 
starburst-generating merger as the "active fraction," eaf- 

Using our merger tree code, we find a relation for the 
active fraction of halos as a function of mass that is nearly 
independent of redshift. Since eaf depends only on the 
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Fig. 1. — The active fraction, eaf, or the fraction of halos 
that have had at least one starburst-generating merger in their 
lifetimes as a function of halo mass, m^, for three different val- 
ues of the galaxy formation suppression threshold mass -A/supp • 
Squares, circles, and triangles show merger-tree simulation results 
for m Bupp = 8, 9.4, and 10, respectively, while dotted, solid, and 
long-dashed lines show the results from equation J9j for the same 
suppression mass values. Enough merger histories we generated so 
that each point represents at least 100 active galaxies. 
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Fig. 2. — The minimum reduced chi-squared (i.e. chi-squared per 
degree-of- freedom) marginalized over Liq as a function of m Bupp . 
Solid, short-dashed, and long-dashed lines show fits to the data 
at z = 6, 7, and 8, respectively. The bottom, middle, and top 
sets of horizontal lines denote the minimum reduced chi-squared 
values required for rejection with 70%, 95%, and 99% confidence, 
respectively, for the number of degrees-of-freedom corresponding 
to the data at each redshift. 



distribution of merger histories, it is also independent 
of e and the other details of our star formation model 
in §2.21 We show, in Figure [U eaf as a function of 
mass calculated from our code for several values of rn supp . 
Each point was generated using enough random merger 
histories to produce at least 100 active galaxies. A good 
fitting formula for caf given by: 



logio eAF("ih,m ; 



1 



suppy 



8 (m h - m sup p) 



2.6 



. (9) 



where and m supp are defined by equation ©. 
Throughout the rest of this paper, we rely on equation © 
to compute eaf for continuous ranges of mh and m supp . 
Our results show that the abundance of halos hosting 
galaxies is suppressed even for halo masses up to an 
order-of-magnitude larger than the suppression thresh- 
old. As we will see in Sj5j the large range of suppression 
masses combines with the distribution of possible lumi- 
nosities for each halo to result in a gentle cut-off at the 
faint end of the LF rather than a sudden drop at a criti- 
cal luminosity. However, our model naturally reproduces 
a high value of 6af for the largest halo masses, consis- 
tent with the rapid evolutio n of the halo mass function 
at z > 6 (|Trenti et al.ll2010[ ). 

5. FITTING THE LUMINOSITY FUNCTION 

We fit our model LF to the late st d ata available at 
z = 6 — 8 from lBouwens et all (|2007ft and lBouwens et all 
(2010c) adopting the same magnitude conventions and 
ignoring, for simplicity, any bright-end upper limits. All 
magnitudes we reference in this paper are rest-frame UV 



absolute magnitudes in the AB system. We have calcu- 
lated LFs at single redshifts for comparison with obser- 
vations, ignoring for the time being the mass-dependent 
distribution of galaxies over the photometric redshift 
range of high-redshift s urveys and its effects on the LF 
(jMunoz fc Loebll2008b| ). At each redshift, we allow two 
fit parameters: L 10 — L c (Mh = 10 10 M Q ) and M supp . 
Liq is directly related to the star formation efficiency e 
in our model for fixed choices of metallicity and IMF. In 
Figure[2]we plot the minimum reduced chi-squared, X?edJ 
(i.e. chi-squared per degree of freedom) values matching 
the observed LF at z = 6, 7, and 8 as a function of M supp . 
Values of Liq have been calculated to minimize %^ cd for 
each value of M supp . 

A minimum in x^ c d a PP e ars at m supp pa 9.5 at z = 6, 
pa 9.4 at z = 7, and « 9.42 at z = 8. At z = 6, all 
combinations of Lin and M SUDD are ruled out at the 70% 

< 8.55 and > 9.7 are 
> 9.8 is ruled out 
8, no constraints 



v 10 c ^ L1, -i J|/J supp 

level. However, values of m supp 
ruled out at the 95% level, while m supp 
with 99% confidence. At z = 7 and z = 
are placed on the minimum value of m supp at the 70% 
level or stronger. However, m supp > 9.7 (9.7), > 9.8 
(9.85), and > 9.9 (9.95) are ruled out with 70%, 95%, 
and 99% confidence, respectively, at z — 7 (8). 

These results clearly indicate that, while the masses of 
halos hosting observed LBGs are typically thought to be 
> 10 10 M Q , lower luminosity galaxies must exist in halos 
smaller than 1O 1O M0, corresponding to a virial temper- 
ature of about 2 x 10 5 K, and very likely in ones at least 
as small as 5 x 10 9 M Q (10 5 K). They also tentatively 
suggest that the minimum mass halo capable of hosting 
galaxies may be around 2.5 x 10 9 M (7 x 10 4 K) with 
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Fig. 3. — A comparison of our best-fit LFs to the data. The 
top, middle, and bottom panels display results for z = 6, 7, and 8, 
respectively. The po i nts and error-bars mark observations from 
IBouwens et al.l {2007, 2010c). Dotted, short-dashed, and long- 
dashed curves are LFs assuming m supp = 8, 9, and 10, respectively, 
with the best-fit value of Lio for each value of m Bupp . Finally, solid 
lines show results with the absolute minimum value of chi-squarc 
at each redshift. The best-fit values of m Bupp are 9.47, 9.4, and 
9.42, for z = 6, 7, and 8, respectively. 



halos less massive than about 3.5 x 10 s M (1.7 x 10 4 K) 
unable to host galaxies with some confidence given the 
data at z = 6. 

Chi-squared is minimized when L\q k, 27.2, 27.4, and 
27.7 at z = 6, 7, and 8, respectively. For our choices of 
metallicity and IMF, these values imply that galaxy for- 
mation is relatively inefficient, with very small fractions 
of galactic gas (0.2%, 0.4%, and 0.5% for each redshift) 
being converted into stars per dynamical time. 

Our best-fit LFs, along with ones assuming m SU pp = 8, 
9, and 10, are shown for each redshift i n Figure [3] The 
data from IBouwens et al.l (|2007l . l2010cD are plotted for 
comparison. The best-fit model deviates qualitatively 
from a Schechter function fit outside the observed mag- 
nitude range. At the bright end, for magnitudes < —21, 
our predicted LF remains much flatter than a Schechter 
fit, which is already beginning to drop exponentially. The 
shallower slope is due to two effects: first, the exponential 
tail of the halo mass function falls off more slowly with 
increasing mass than a Schechter function with luminos- 
ity proportional to mass, and second, the large spread 
in the luminosity permitted for each halo mass allows 
abundant, smaller halos emitting at higher than average 
luminosity to bolster the population of bright galaxies. 

On the other hand, the suppression of star formation 
drastically reduces the abundance of galaxies at mag- 
nitudes fainter than currently observable compared with 
expectations from a simple extrapolation of the Schechter 
function. The result is a flatter LF slope in the ob- 
served region for increasing m supp . Figure [3] clearly il- 
lustrates the disparate predictions for the abundance of 



faint galaxies between different fiducial values of M supp . 
Additional data at only about a magnitude fainter than 
the current observational threshold will greatly improve 
our ability to constrain the minimum halo mass capable 
of forming galaxies. For reference, observations down to 
a magnitude of -16.8 at z — 7 will require a sensitivity 
of about 1.5 nJy, close to what is expected with JWST. 
However, while the lcr errors in the current data were cal- 
culated based on the shot noise and cosmic variance from 
an amalgam of observations from several different fields, 
we conservatively estimate the la error on the abundance 
at this magnitude in a single 2' x 2' pointing of NIRCam 
on JWST to be about 50% (IMunoz et al.ll201(l . 

We have explicitly ignored the influence of quiescent 
star formation in our model. Such a mechanism in small 
halos may reduce the effects that we describe of a galaxy 
suppression mass on the faint end of the LF making it 
more difficult to pr obe such phy s ics wi th future surveys. 
However, work by lLacev et al.l (|2010l ) has shown that 
merger-driven starbursts do dominate the UV LF down 
to at least magnitude -17 at z—6 and to -15 by z=10, 
albeit with a very different IMF. Thus, while more com- 
plicated models may be required for high-precision mea- 
surements of -M supp even with a complete LF, we are 
confident that deeper surveys of the not-to-distant fu- 
ture will help illuminate some of the physics of galaxy 
suppression. 

6. STAR FORMATION RATE 

The SFR of high-redshift galaxies is important for un- 
derstanding the sta r formation history of the Universe 
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1998) and the ionization state of the IGM 

1999) . Its estimation relies on a propor- 



tionality between UV luminosity and SFR based on two 
assumptions: an exponential burst of star formation has 
a timescale, t, that is longer than 1 Gyr, and the stellar 
populatio n is observed after one exponential time scale 
has past (jMadau et al.|[l998l ) . However, if the age of the 
universe is shorter than 1 Gyr, at least one of these as- 
sumptions must be violated. 

For the best-fit star formation efficiencies we found 
from the data at z = 6, 7, and 8, the typical exponen- 
tial starburst timescale given by equation §5§ is of order 
t 10 Gyr, an order-of-magnitude or more longer than 
the age of the universe. Equation gives the typical 
SFR to be of order lMQ/yr in a burst with 10 10 M Q 
worth of gas remaining; if the amount of initial gas in a 
halo as a fraction of the total halo mass is S!b/^ m ~ 0.16, 
this corresponds to the initial SFR in a halo of about 
6.25 x 10 10 M Q . 

Figured] shows the evolution of the luminosity at 1500 
-^1500 1 an d SFR, M*, with time and their relation- 
ship calculated for exponential bu rsts using a simple ste l- 
lar population from Starburst99 (jLeitherer et al.l fl999 h 
Solid lines represent typical bursts in high-redshift galax- 
ies with the initial SFR set at 1 M / yr and the exponen- 
tial time scale r = 10 Gyr. For t > r, both M* and L1500 
decrease exponentially over time with timescale r so that 
-£1500 is proportional to M+. This is because the exponen- 
tial timescale is much longer than the lifetime of the stars 
that dominate the UV luminosity. The amplitude of the 
relation is set by the IMF and metallicity of the stellar 
population; for the choices described in ^2.2i we find ap- 
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10 6 10 7 10 6 10 9 10 10 10 11 

burst age / yrs 

Fig. 4. — The UV luminosity and SFR evolution of exponen- 
tial bursts with (r/Gyr, M+(t = Q)M -1 yr) = (10,1), (0.1,3), 
and (0.01,10) denoted by solid, long-dashed, and short-dashed 
curves, respectively. The top panel tracks the bursts in SFR- 
^1500 space. Here, the upper and lower dotted lines show a pro- 
portional relationship between SFR and L1500 with constants of 
8 X 10 27 erg/s/Hz/(M Q /yr) and 2 X 10 28 erg/s/Hz/(M /yr), re- 
spectively. The middle panel shows the burst lightcurves with right 
and left vertical lines denoting the r = 10 Gyr and the age of the 
universe at z = 7, respectively. The lightcurve for an instantaneous 
burst producing 10 6 Mg worth of stars is given by the dotted curve. 
The horizontal line marks the observable threshold at a magnitude 
of -18. The bottom panel plots the evolution of the SFR with time, 
while the dotted curve here shows the SFR expected from by the 
relationship between SFR and L1500 given the burst luminosity as 
a function of time. 



proximately £1500 — 2 x 10 28 (M*/M Q yr) erg/s/Hz with 
a proportionality constant a factor of 2.5 different than 
the 8 x 10 27 erg/s/Hz/(M Q /yr). However, before t = t, 
the luminosity is still rising with increasing time, while 
the SFR remains essentially unchanged. Since the age of 
the universe is much less than r, all bursts are observed 
in this phase before the £1500 — M* proportionality has 
stabilized. 

Thus, the SFR will typically be somewhat higher than 
that inferred from the £1500 — proportionality. The 
ratio between the true and expected SFRs will depend 
on how close the burst is to its maximum luminosity, the 
point where the expected SFR is approximately equal 
to its initial value. The burst lightcurve is relatively 
flat near its maximum value over the time approximately 
10 7 — 10 9 yrs after it begins. If the burst is more than 
10 7 yrs old at the time of observation, the difference be- 
tween the true and expected SFRs will not be very sig- 
nificant. A burst observed at z = 7 is 10 7 yrs old if it 
started at z w 7.07. 

For completion, we also show in Figure [4] the evo- 
lution of L1500 and M* for bursts with r = 0.1 and 
0.01 Gyr, les s than the 1 Gyr minimum considered by 
iMadau et al.l (|1998| ). These timescales are achieved at 
z = 7 if for combinations of the star formation efficiency 




-22 -20 -18 -16 -14 



Fig. 5. — The SFR vs. rest-frame UV magnitude of simulated 
halos. Magenta, blue, green, and red points denote halos of mass 
10 10 , 10 10 ' 5 , 10 11 , and lO n - 5 M , respectively. The solid, black 
line marks L1500 = 2 X 10 2s (M*/M yr) erg/s/Hz. 

and the disk spin parameter such than A = 0.17 and 
0.017, respectively. The luminosity in each case begins 
to decline before reaching the maximum it would have 
achieved had r been longer. The fall-off in luminosity is 
only slightly slower than exponential for r = 0.1 Gyr so 
that the SFR and UV luminosity reach a nearly propor- 
tional relationship after t — t, albeit with a coefficient 
slightly higher than that seen for higher r. However, the 
decline in luminosity is more power-law than exponential 
for r = 0.01 Gyr leading to a very non-linear relationship 
between SFR and luminosity after t = r. In both cases, 
the SFR is much less than expected for a given lumi- 
nosity. This is because the timescale r is not so much 
longer than the lifetimes of the stars that dominate the 
UV luminosity. Luminosity from stars produced at t = r, 
for example, contribute significantly to the luminosity at 
0.1 r, whereas the contribution would be completely neg- 
ligible for much larger r. Consequently, the luminosity 
for a given instantaneous SFR can be much higher than 
expected. 

Using our merger tree and star formation code, we cal- 
culate the instantaneous M* at the time of observation 
for ea ch of our modeled g alaxies and test the accuracy 
of the IMadau et alj () 19981 ) proportionality at z = 7 over 
a wide range of halo masses and for a full distribution 
of spin parameters and merger histories. Figure [5] shows 
the relationship between instantaneous SFR, M*, and 
UV luminosity, L1500J for galaxies in halos at z — 7 with 
masses of 10 10 , 10 10 5 , 10 11 , and 1O 115 M . We have set 
Wsupp — 9.4. Each point represents a single halo, and 
where more than one ongoing starburst is present, we 
have simply added the contributing SFRs. The discrete 
" lines" of points above the main body for the smaller halo 
masses is a resolution effect of our code; while their exact 
positions should not be taken as precise, such points do 
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represent a population of halos lying above the standard 

^1500 — M* relation. 

The majority of points in Figure [5] do show a rough 
proportionality between L1500 and M+. However, the 
proportionality constant is slightly higher than the 2 x 
10 28 erg/s/Hz/(M Q /yr) value for bursts with ages longer 
than their exponential time scale with the difference de- 
pending on halo mass. Lower mass halos tend to be pop- 
ulated by younger bursts that are further from reaching 
their maximum luminosity than higher mass halos. If we 
constrain £1500 <x M*, we find a proportionality constant 
of 1.7 x lO 28 erg/s/Hz/(M /yr) for 1O 1O M halos which 
estimates SFRs to be about 15% higher than the con- 
stant for older bursts. Given the typical uncertainties in 
measuring the total star formation rate at high redshift 
- sample completeness, cosmic variance, uncertain IMF 
and metallicity, etc. - an additional ~ 20% error will not 
significantly affect current estimates. 

7. IONIZATION STATE OF THE IGM 

After cosmic reionization, the ionization state of the 
IGM depends on the balance between the recombination 
rate and the production rate of ionizing photons. On its 
own, the formation of stars in galaxies can maintain the 
ionization of the IGM through its production of ioniz- 
ing photons if the star formation ra te density (SFRD) i s 
higher than a critical value given bv lMadau etaLl (|1999() : 



2xW- 3 f- 1 c C 



1 + z 
10 



M /yr/Mpc 3 , (10) 



where / csc is the fraction of ionizing photons produced 
in galaxies that escape into the IGM, and C is the IGM 
clumping factor. 

Using the standard £1500 to conversion, recent ob- 
servational studies have estimated the currently observ- 
able SFRD to be just enough to keep the universe ionized 
if f~c C = 1. However, much of the star formation be- 
low the observable threshold is not included. In Figure 
[BJ we compare our calculations for the total SFRD at 
z = 6, 7, and 8 assuming the best-fit values of L w and 
M aupp at each redshift to the observed estimates and to 
equation (flU)) . We show results for L 150 o to ratios 
of both 8 x 10 27 erg/s/Hz/(M Q /yr) (the typically used 
value) and 2 x 10 28 erg/s/Hz/(M /yr) (consistent with 
our choices of IMF and metallicity). We also show the 
factor by which the total values of the SFRD and UV lu - 
min osity exceed those obse rved by lBouwens et al.l (|2007f ) 
and lBouwens et all (|2010cfl . Factors less than unity indi- 
cate that the observed points are higher than average in 
the universe at that redshift due to Poisson fluctuations 
or cosmic variance so that the observed SFRD is higher 
than the average over the whole population. 

Our results show that the ability of the total galaxy 
population to account for the UV background required 
to keep the IGM ionized depends critically on the value 
of m supp . For m supp = 8, the total SFRD or UV lumi- 
nosity is about 3-9 times the observed values with more 
star formation and luminosity missing at higher redshift. 
However, at the best-fit values of m supp = 9.5 at z — 6 
and 9.4 at z = 7 and 8, the galaxy population produces 
no more SFRD than observed (and somewhat less for 
z = 6). Assuming /~* C = 1, the total SFRD for all pa- 




FlG. 6. — The top panel shows the SFRD produced by the to- 
tal galaxy population at z = 6, 7, and 8. Circles denote re- 
sults using best-fit values of M sup p and L10 at each redshift, 
while triangles assume Af sup p = 10 8 Mq and the corresponding 
best-fit values of Lio- Filled (empty) points use Lisoo/M* = 
2 X 10 28 (8 X 10 27 )erg/s/Hz/(Mp)/yr ) Square poin t s wit h er- 
ror bars denot e observ ed values from Bouwcns et al. (2003) and 
IBouwens et~a l. (2010c). The minimum SFRD required to keep 
the IGM ionized as given by Eq. I|10ll for f^sc C = 15 and 1 are 
shown by the upper and lower solid lines, respectively. The bot- 
to m panel shows the rati o of the total UV luminosit y or SFRD to 
the Bouwcns et al. (2007) and Bouwcns ct al. (2010c) observations 
as a function of M supp . The solid, short-dashed, and long-dashed 
lines denote z = 6, 7, and 8, respectively. 

rameters and redshifts considered meet the requirement 
for maintaining the ionization of the IGM. However, if 
/esc C =15 (e.g. / osc = 0.2 and C = 3), the galaxy popu- 
lation can keep the IGM at z = 8 ion ized only for a choice 
of IMF and metallicity that gives the lMadau et al.1 ()1998f ) 

ratio of Li 500 to M* ratio of 3 x 10 27 erg/s/Hz/(M /yr) 
and if m supp ~ 8, lower than our best-fit value. Finally, 
since the amount of star formation below the observable 
limit increases with redshift, we find that the evolution 
of the total SFRD with redshift is much flatter than that 
observed. 

8. CONCLUSIONS 

In this paper, we combine a standard merger tree algo- 
rithm with a simple star formation prescription designed 
to encapsulate the main physical processes relevant at 
z > 6. Our model both accounts for a range of possi- 
ble galaxy luminosities for each halo mass and includes a 
sharp galaxy formation cut-off in halo mass below Af supp . 

• We confirm that the luminosity distribution func- 
tion for halos of a given mass is a roughly log- 
normal distribution with a variance of ~ 1.5 magni- 
tudes and a proportional relationship between the 
mean luminosity and halo mass (see 



At a fixed halo mass of 10 M , 
luminosities are logio(£io erg" 1 s Hz) = 



the mean log- 
27.2, 27.4, and 
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27.7 at z = 6, 7, and 8, respectively, suggesting that 
for a fixed halo mass, galaxies are brighter on average 
at higher redshift, consistent with results from Schechter 
fits. However, while the exponential tail of the high- 
redshift halo mass function is less sharp than that of a 
Schechter function, the range of possible luminosities for 
a fixed halo mass further slows the fall-off of the predicted 
galaxy LF at the bright end. While still being consistent 
with the data, our shallower LF anticipates the discov- 
ery of a larger population of very bright galaxies at z = 7 
and 8 as survey fields increase in size. 

• We also show that an active fraction of halos that 
approaches unity with increasing halo mass can be 
naturally explained by a suppression halo mass for 
galaxy formation combined with the variety of pos- 
sible merger histories (see $4]). 

This active fraction is well-approximated by the for- 
mula given in Eq. ([9]). One can easily use this expression, 
along with our log-normal distributions of UV luminos- 
ity for each halo mass to calculate the galaxy LF from 
the halo mass function. The resulting LF does not have 
a sharp cutoff at the faint end but rather turns over gen- 
tly. Thus, we predict that as long as future observations 
show a LF that increase with ever decreasing luminosity, 
the surveyed region will never be volume-complete. 

• The current data suggests that the minimum mass 
halo capable of hosting galaxies may be around 
2.5 x 1O 9 M , corresponding to a virial tempera- 
ture of 7 x 10 4 K (see 

We find a strong upper limit of M supp < 6 x 10 9 M 
(10 5 K) with at least 95% confidence. However, lower 
limits from the current data are quite weak with halos 
less massive than about 3.5 x 10 8 M (1.7 x 10 4 K) unable 
to host galaxies with some confidence given the data at 
z = 6. 

• We find a best-fit star formation efficiency at high 
redshift of approximately 0.2-0.5% per dynamical 
time, implying a starburst exponential time scale 
much longer than the age of the universe. 

A more top-heavy IMF would have required even less 
efficient star formation, corresponding to even longer 
burst time scales, to produce the same observed lumi- 
nosities. However, the long burst time scale does not 
create lightbulb-like galaxies that, once switched on, are 
always emitting with the same luminosity. Instead con- 
tinued merger activity disrupts old bursts and replaces 
them with new ones based on the particular history of 
the host halo. 

• We show that the proportionality of -L1500 to M* is 
usually an adequate approximation (see $6j. 



While the iMadau et al.l (|T998) proportionality relies 
on long-lived bursts in the tail of their exponentially de- 
creasing rate of star formation, most bursts at z = 7 are 
emitting near their maximum luminosity where the track 
of SFR vs. UV luminosity begins to join the proportional 
relationship. Despite their young ages compared to their 
exponential time scale, this is because the bursts are typ- 
ically older than 10 7 yrs at the time of observation, old 
enough that the massive stars providing the bulk of the 
UV luminosity are beginning to die out as fast as new 
ones are added. Although the lowest mass halos may 
host very young bursts that have somewhat higher SFRs 
than expected for their luminosities, using a standard 

proportionality of £1500 to M* adds additional errors of 
only tens of percent. However, some care must be taken 
in selecting a constant of proportionality consistent with 
spe cific choices of metallicity and IMF rather than using 
the IMadau et al.1 (|1998l) value indiscriminately. Addi- 
tionally, since bursts are likely to remain close to their 
maximum SFRs and luminosities for most of their life- 
times, ongoing accretion between major mergers is less 
likely to be important. 

• When extrapolated down to faint luminosities be- 
low the current observable threshold, the total 
SFRD of the galaxy population will only at most 3- 
9 times higher than what has already been observed 
even if the minimum halo mass forming galaxies is 
as low as 1O 8 M (see 

The gentle turnover at the faint end of the LF, even 
given a sharp cutoff in the halo mass capable of producing 
galaxies, results in less star formation below the observ- 
able limit than if the LF dropped sharply at the mean 
luminosity corresponding to the same halo mass. While 
the total galaxy population with m supp = 8 may be able 
to keep the IGM ionized given /~* C ~ 15, for our best- 
fit value of TO sup p « 9.4, no significant star formation lies 
below a rest-frame UV magnitude of -18. In such a case, 
galaxies may only be responsible for maintaining the ion- 
ization state of the IGM if /~* C ~ 1 . Interestingly, we 
also find that, since the amount of missing star forma- 
tion increases with redshift, the redshift evolution of the 
total star formation history of the universe is flatter than 
observed. 

Although the current data from LBG drop-outs does 
not place a strong lower-limit on the minimum halo mass 
required to host galaxies at redshifts > 6, we have shown 
that JWST and other future deep surveys will provide 
much tighter constraints. These results will not only 
shed light on the contribution of galaxies to the UV back- 
ground that keeps the IGM ionized but also hint at the 
feedback physics that limits galaxy formation. 

9. ACKNOWLEDGEMENTS 

We thank Steve Furlanetto for useful discussions. This 
work was supported in part by NSF grant AST-0907890 
and NASA grants NNX08AL43G and NNA09DB30A 
(for A.L.). 



REFERENCES 



Barkana, R., & Loeb, A. 2000a, ApJ, 531, 613 
— . 2000b, ApJ, 539, 20 



Baugh. C. M., Lacey, C. C, Frenk, C. S., Granato, G. L., Silva, 
L., Bressan, A., Benson, A. J., & Cole, S. 2005, MNRAS, 356, 
1191 



Minimum Galaxy Mass at High-z 



9 



Bouwons, R. J., Illingworth, G. D., Franx, M., & Ford, H. 2007, 

ApJ, 670, 928 
Bouwens, R. J., et al. 2008, ApJ, 686, 230 
Bouwcns, R. J., Illingworth, G. D., Oesch, P. A., Trenti, M., 

Stiavclli, M., Carollo, C. M., Franx, M., van Dokkum, P. G., 

Labbe, I., & Magee, D. 2010a, ApJ, 708, L69 
Bouwens, R. J., Illingworth, G. D., Labbc, I., Oesch, P. A., 

Carollo, M., Trenti, M., van Dokkum, P. G., Franx, M., 

Stiavelli, M., Gonzalez, V., & Magee, D. 2010b, ArXiv e-prints 
Bouwens, R. J., Illingworth, G. D., Oesch, P. A., Labbe, I., 

Trenti, M., van Dokkum, P., Franx, M., Stiavelli, M., Carollo, 

C. M., Magee, D., & Gonzalez, V. 2010c, ArXiv e-prints 
Cole, S., Lacey, C. G., Baugh, C. M., & Frenk, C. S. 2000, 

MNRAS, 319, 168 
Cooray, A., & Milosavljevic, M. 2005a, ApJ, 627, L85 
— . 2005b, ApJ, 627, L89 

Cooray, A., & Ouchi, M. 2006, MNRAS, 369, 1869 

Dunklcy, J., Komatsu, E., Nolta, M. R., Spergel, D. N., Larson, 

D. , Hinshaw, G., Page, L., Bennett, C. L., Gold, B., Jarosik, 
N., Weiland, J. L., Halpern, M., Hill, R. S., Kogut, A., Limon, 
M., Meyer, S. S., Tucker, G. S., Wollack, E., & Wright, E. L. 
2009, ApJS, 180, 306 

Fakhouri, O., Ma, O, & Boylan-Kolchin, M. 2010, MNRAS, 406, 
2267 

Finlator, K., Oppenheimer, B. D., & Dave, R. 2010, ArXiv 
e-prints 

Gonzalez, V., Labbe, I., Bouwens, R. J., Illingworth, G., Franx, 

M., Kriek, M., & Brammer, G. B. 2010, ApJ, 713, 115 
Kennicutt, Jr., R. C. 1998, ApJ, 498, 541 

Komatsu, E., Dunkley, J., Nolta, M. R., Bennett, C. L., Gold, B., 
Hinshaw, G., Jarosik, N., Larson, D., Limon, M., Page, L., 
Spergel, D. N, Halpern, M., Hill, R. S., Kogut, A., Meyer, 
S. S., Tucker, G. S., Weiland, J. L., Wollack, E., & Wright, 

E. L. 2009, ApJS, 180, 330 

Labbc, I., Gonzalez, V., Bouwens, R. J., Illingworth, G. D., Oesch, 
P. A., van Dokkum, P. G., Carollo, C. M., Franx, M., Stiavclli, 
M., Trenti, M., Magee, D., & Kriek, M. 2010a, ApJ, 708, L26 

Labbc, I., Gonzalez, V., Bouwens, R. J., Illingworth, G. D., Franx, 
M., Trenti, M., Oesch, P. A., van Dokkum, P. G., Stiavelli, M., 
Carollo, C. M., Kriek, M., & Magee, D. 2010b, ApJ, 716, L103 



Lacey, C. G., Baugh, C. M., Frenk, C. S., Benson, A. J., & . 

2010, ArXiv c-prints 
Leitherer, C, Schaerer, D., Goldader, J. D., Gonzalez Dclgado, 

R. M., Robert, C, Kune, D. F., de Mello, D. F., Devost, D., & 

Heckman, T. M. 1999, ApJS, 123, 3 
Madau, P., Haardt, F., & Rees, M. J. 1999, ApJ, 514, 648 
Madau, P., Pozzetti, L., & Dickinson, M. 1998, ApJ, 498, 106 
Mo, H. J., Mao, S., & White, S. D. M. 1998, MNRAS, 295, 319 
Munoz, J. A., & Locb, A. 2008a, MNRAS, 385, 2175 
— . 2008b, MNRAS, 386, 2323 

Munoz, J. A., Trac, H., & Loeb, A. 2010, MNRAS, 405, 2001 
Oesch, P. A., Bouwens, R. J., Carollo, C. M., Illingworth, G. D., 

Trenti, M., Stiavelli, M., Magee, D., Labbe, I., & Franx, M. 

2010, ApJ, 709, L21 
Ovcrzicr, R. A., Guo, Q., Kauffmann, G., De Lucia, G., Bouwens, 

R., & Lemson, G. 2009, MNRAS, 394, 577 
Raicevic, M., Theuns, T., & Lacey, C. 2010, ArXiv e-prints 
Robertson, B. E. 2010a, ApJ, 713, 1266 
— . 2010b, ApJ, 716, L229 

Salvaterra, R., & Ferrara, A. 2006, MNRAS, 367, Lll 
Salvaterra, R., Ferrara, A., & Dayal, P. 2010, ArXiv e-prints 
Somerville, R. S., & Kolatt, T. S. 1999, MNRAS, 305, 1 
Stark, D. P., Loeb, A., & Ellis, R. S. 2007, ApJ, 668, 627 
Trenti, M., & Stiavelli, M. 2008, ApJ, 676, 767 
Trenti, M., Stiavelli, M., Bouwens, R. J., Oesch, P., Shull, J. M., 

Illingworth, G. D., Bradley, L. D., & Carollo, C. M. 2010, ApJ, 

714, L202 

Vale, A., & Ostriker, J. P. 2004, MNRAS, 353, 189 
— . 2006, MNRAS, 371, 1173 
— . 2008, MNRAS, 383, 355 

Wetzel, A. R., Cohn, J. D., & White, M. 2009, MNRAS, 395, 1376 
Wetzel, A. R., & White, M. 2010, MNRAS, 403, 1072 
Wyithe, J. S. B., & Loeb, A. 2006, Nature, 441, 322 



